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Abstract 

Civil  defense  shelters  are  often  constructed  beneath  the  ground  to  provide  protection  against  blast  loads.  Concrete  is  widely  used  as 
the  material  for  the  defense  layer  of  the  shelters  and  artificial  cavities  are  often  embedded  in  the  concrete  to  ensure  adequate  attenuation 
of  the  air-blast  waves.  This  paper  employs  the  Johnson-Holmquist-Concrete  model,  an  elastoplastic  damage  model,  to  investigate  the 
effects  of  artificial  cavities  on  the  attenuation  of  air-blast  waves  in  the  concrete  defense  layer.  The  results  reveal  that  the  peak  axial  stress 
beneath  a  cavity  is  largely  reduced  and  the  elastoplastic  and  damage  properties  of  the  concrete  play  significant  roles  on  the  wave  atten¬ 
uation.  Two  empirical  formulae  are  then  proposed  to  relate  the  decay  factor  of  peak  axial  stress  to  the  dimensions  and  relative  positions 
of  rectangular  and  circular  cavities.  Finally,  the  effects  of  multiple  cavities  on  the  attenuation  of  air-blast  waves  are  investigated. 

©  2006  Elsevier  Ltd.  All  rights  reserved. 
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1.  Introduction 

The  propagation  and  diffraction  of  air-blast  induced 
stress-waves  are  important  considerations  in  anti-blast 
design.  The  development  of  modern  military  technology 
greatly  enhances  the  destructive  power  and  hit  rate  of 
new  weapons,  thus  posing  great  challenges  to  defense  engi¬ 
neering  nowadays.  Civil  defense  shelters  are  widely  con¬ 
structed  beneath  the  ground  to  provide  protection  against 
blast  loads  [1,2].  Blast  waves  due  to  far-held  detonation 
evolve  into  propagation  of  stress-waves  in  underground 
media.  Such  waves  will  diffract  when  hitting  obstacles,  such 
as  cavity,  crevice  or  other  media.  The  wave  stress  intensity 
may  be  greatly  reduced  beneath  these  obstacles.  This 
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reduction  of  stress  from  the  peak  value  is  often  termed  as 
attenuation  or  insulation  of  stress-waves  [3,4]. 

The  diffraction  and  attenuation  of  stress-waves  are 
important  aspects  in  the  design  of  civil  defense  shelters 
[1,4-6].  To  design  a  civil  defense  structure,  the  reduction 
of  stress-waves  needs  to  be  examined.  In  addition,  new 
stress-wave  migrating  and  weakening  systems  are  expected 
[7].  In  this  paper,  investigations  will  be  focused  on  the 
reduction  of  stress-waves  through  an  intelligent  defense 
layer  embedded  with  artificial  cavities. 

Fig.  1(a)  shows  that  an  intelligent  defense  layer  usually 
consists  of  three  layers  [5,8]:  an  earth  cover  layer,  a  pro¬ 
tection  layer  and  a  support  layer.  The  protection  layer 
often  has  two  sub-layers:  a  projectile  shelter  layer  and  a 
stress  distribution  layer.  The  projectile  shelter  layer  pro¬ 
vides  resistance  against  penetration,  and  several  design 
schemes  on  the  topic  have  been  proposed  so  far  [8-10]. 
The  stress  distribution  layer  is  employed  to  redistribute 
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2.  Computational  problem  for  a  typical  intelligent  layer 
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Fig.  1 .  Typical  defense  shelter  and  simulation  of  stress  distribution  layer, 
(a)  Typical  defense  shelter,  (b)  Computational  domain  for  stress  distribu¬ 
tion  layer. 


or  diffuse  the  blast  loading.  Till  to  date,  relatively  few 
studies  have  been  carried  out  to  examine  the  stress  redis¬ 
tribution  in  this  layer.  For  example,  artificial  cavities  are 
constructed  in  the  stress  distribution  layer  to  dissipate 
the  energy  of  stress-waves.  However,  the  function  and 
design  methodology  of  these  cavities  have  not  been  thor¬ 
oughly  examined  so  far  [11-16].  Furthermore,  concrete  is 
commonly  employed  as  the  material  of  these  intelligent 
defense  layers  in  protective  facilities.  Upon  blast  loading, 
the  concrete  medium  usually  experiences  elastoplastic 
deformation  and  even  damage  and  this  may  affect  the 
attenuation  of  the  stress- waves.  However,  the  effects  of 
elastoplastic  and  damage  properties  of  concrete  on  the 
attenuation  of  stress-waves  are  not  well  understood.  In 
addition,  very  little  is  known  on  the  screen  effect  of 
cavities. 

In  this  paper,  numerical  analyses  have  been  carried  out 
to  investigate  the  effect  of  an  embedded  cavity  on  the  atten¬ 
uation  of  stress-waves  in  concrete  medium  using  the  John- 
son-Holmquist-Concrete  (J-H-C)  model  [17,18].  In 
particular,  the  influence  of  cavity  on  the  decay  factors  of 
peak  axial  stress,  oy ,  is  examined  in  detail  using  rectangular 
and  circular  shaped  cavities.  Based  on  the  numerical 
results,  empirical  formulae  are  proposed  to  relate  the  decay 
factor  for  peak  stress  to  the  cavity  parameters.  Finally,  the 
effects  of  multiple  cavities  are  studied  and  compared  with 
that  of  single  cavity. 


To  evaluate  the  attenuation  extent  of  a  stress-wave  due 
to  a  cavity  in  an  intelligent  defense  layer,  a  dimensionless 
variable,  'decay  factor’  DF,  is  defined  as  follows: 


DF  = 


do  -  a 


(Jo 


(i) 


where  a0  denotes  the  peak  value  of  stress  component,  such 
as  ox ,  Oy,  and  hydrostatic  pressure  P  at  a  specified  position 
without  cavity  in  the  defense  layer,  and  o  is  the  peak  value 
of  the  same  component  at  the  same  position  with  a  cavity 
in  the  defense  layer.  Eq.  (1)  reveals  that  a  larger  DF  would 
imply  a  superior  attenuation  effect. 


2.2.  Computational  domain  with  circular  cavity 

A  typical  concrete  defense  layer  embedded  with  a  circu¬ 
lar  cavity  is  shown  in  Fig.  1(b).  This  is  a  two-dimensional 
model  whose  computational  domain  is  taken  as 
600  cm  x  300  cm  (width  x  height).  The  upper  apex  of  the 
cavity  is  at  a  depth  e  beneath  the  surface  and  the  circular 
cavity  has  a  radius  R.  Thus,  the  coordinates  of  the  center 
is  (0,  e  +  R).  Other  dimensions  are  labeled  in  Fig.  1(b). 
Introducing  the  term  relative  distance  y  as  the  distance 
between  point  A  and  the  lower  apex  of  the  cavity,  the 
dimension  Y  follows  the  following  relationship: 

Y  =  y  +  2R  +  e  (2) 

Only  half  of  the  domain  is  selected  in  the  computations 
owing  to  symmetry.  As  an  example,  the  right  half  domain 
shown  in  Fig.  1(b)  has  the  following  boundary  conditions: 
pressure  boundary  on  the  top,  transmission  boundaries  (no 
reflection)  on  the  right  and  lower  sides,  and  symmetric 
boundary  on  the  left  side. 

2.3.  Air-blast  loading 

A  spherical  front  of  blast  waves  are  reported  to  be 
developed  during  an  explosion  in  the  air  [19,20].  The  ini¬ 
tial  high  pressure  decreases  as  the  wave  propagates 
through  the  air.  As  the  blast  wave  front  is  much  larger 
than  the  dimensions  of  the  defense  layer,  the  air-blast 
loading  or  overpressure  can  be  assumed  to  be  uniform 
over  the  entire  upper  surface  of  the  defense  layer.  In 
addition,  because  the  overpressure  decays  exponentially, 
it  can  be  approximated  by  an  equivalent  triangular  pres¬ 
sure  pulse  as  shown  in  Fig.  2  [21,22].  The  intensity  and 
duration  of  this  overpressure  depend  on  the  explosive 
charge  weight  and  the  stand-off  distance  which  usually 
follows  the  cubic  laws  [1,21-23].  In  the  present  study, 
the  peak  value  Qmax  is  taken  as  120  MPa  and  occurs 
at  time  tmax  =  0  (is  (micro  second).  This  peak  overpres¬ 
sure  then  decreases  linearly  to  zero  at  time  tc  —  100.0  (is. 
Such  an  overpressure  corresponds  to  the  charge  weight 
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o  tc 


Fig.  2.  Triangular  impulse  loading  for  explosion. 

of  885  kg  (TNT  explosive)  and  stand-off  distance  of  5  m 
in  the  air. 

2.4.  Numerical  tools  for  analysis 

The  LS-DYNA  software  from  the  Livermore  Software 
Technology  Corporation  is  used  in  the  present  study  to 
carry  out  all  computations.  This  software  can  handle  either 
two-  or  three-dimensional  problems  for  high-speed  impact 
and  explosion  [24,25]. 

3.  Constitutive  law  of  concrete  medium 

Concrete  is  typically  nonhomogeneous  in  composition 
and  may  possess  flaws,  such  as  cracks  and  voids.  Strictly 
speaking,  the  linear-elastic  model  (L-E  model)  can  only 
describe  its  initial  deformation  [17,26].  Under  air-blast 
loading,  the  flaws  may  cause  elastoplastic  and  damage 
deformations  in  the  concrete  medium.  The  elastoplastic 
damage  process  is  irreversible  and  dissipates  the  stress- 
wave  energy.  This  paper  adopts  the  J-H-C  model 
[17,18,23],  which  has  been  specifically  developed  for  con¬ 
crete  subject  to  high  pressures,  large  strains,  and  high  strain 
rates.  As  seen  from  Fig.  3,  this  model  has  three  compo¬ 
nents:  strength  criterion  with  strain  rate  effect,  damage  evo¬ 
lution  with  plastic  deformation,  and  nonlinear  hydrostatic 
compression. 

3.1.  Strength  with  strain  rate  effect 

If  the  equivalent  stress  o  is  normalized  by  the  quasi-sta¬ 
tic  uniaxial  compressive  strength  ff  a  normalized  equiva¬ 
lent  stress  c*  is  defined  as 


a*  is  a  function  of  pressure  P,  strain  rate  s  and  damage 
parameter  D 

a*  =  [A(  1  -D)  +B-P*n}[1  +  CMn(e*)]  (4) 

where  P*  =  j  denotes  the  normalized  pressure,  and  s*  =  4 
is  the  dimensionless  strain  rate.  In  addition,  A  is  the  nor¬ 
malized  cohesion,  B  is  the  coefficient  for  the  normalized 


b 


c 

Fig.  3.  Elastoplastic  damage  constitutive  model  of  concrete  medium, 
(a)  Constitutive  model,  (b)  Damage,  (c)  Hydrostatic  compression  pressure 
versus  volumetric  strain. 

pressure  hardening,  C  denotes  the  coefficient  of  strain  rate, 
and  N  is  an  exponent  for  the  pressure  hardening.  The  dam¬ 
age  parameter  D  will  be  examined  next. 
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3.2.  Damage  parameter 

In  the  present  model,  the  damage  parameter  D  is  defined 
by  the  accumulation  of  both  generalized  deviatoric  plastic 
strain  sp  and  plastic  volumetric  strain  /ip.  It  is  expressed  as 

D  =  y  Agp+A\  (5) 

^  Di(P*  +  T*)Dl 

where  A  denotes  the  increment  and  denotes  their  sum¬ 
mation  along  the  deformation  history;  T*  =  L  is  the  nor¬ 
malized  hydrostatic  pressure;  and  D\  and  D2  are  the 
material  constants.  It  is  noted  that  the  damage  parameter 
D  varies  within  the  range  of  0-1.  When  D  =  0,  no  damage 
is  observed;  and  when  D=  1,  the  concrete  is  completely 
damaged  and  cannot  resist  any  shear  stress. 


3.3.  Hydrostatic  compression 


Since  concrete  is  formed  by  granular  materials,  its 
mechanical  behavior  depends  on  the  hydrostatic  pressure. 
The  hydrostatic  pressure-volume  relationship  is  presented 
in  Fig.  3(c).  The  compressive  deformation  of  concrete  expe¬ 
riences  the  following  three  stages: 


(a)  Elastic  stage  (0  ^  ^  /Wh) 

P  =  Klfi,  n  =  ~,  *i=  —  (6) 

P  0  /Lrush 

where  p ,  p0  are  the  current  and  initial  densities, 
respectively;  p crush  is  the  crushing  volumetric  strain; 
^cmsh  the  crushing  pressure;  and  K\  the  bulk 
modulus. 

(b)  Transition  stage  (/Wsh  <  p  ^  ^iock)  The  piock  is  the 
locking  volumetric  strain.  The  bubbles  in  concrete 
are  broken  in  this  stage,  causing  damage  of  the  con¬ 
crete  structure,  but  the  integrity  of  concrete  medium 
is  still  maintained. 

(c)  Compact  stage  (p  >  p\ock)  The  concrete  is  now  fully 
dense  material  without  air  and  voids.  The  pressure 
P  for  this  fully  dense  material  is  formulated  as 

P  =  Kln  +  K2Ji2  +K3Ji3  (7) 


where  K2 ,  K3  are  material  constants.  The  modi¬ 
fied  volumetric  strain  p  is  defined  as: 


-  P  Plock 

P  i  , 

1  '  Mlock 


(8) 


4.  Numerical  investigations 

4.1.  Circular  cavity 

4.1.1.  Effect  of  cavity  radius  R 

Five  radii  are  chosen  in  the  computations:  0.0,  0.15, 
0.30,  0.45  and  0.60  m,  respectively.  The  parameter  e  is  fixed 
as  0.30  m.  The  model  parameters  for  concrete  are  given  in 


Table  1 

Model  parameters  for  concretea 

Density 

Po  (kg/m3) 

2440 

Elastic  constants 

G  (GPa) 

14.86 

V 

0.20 

K\  (GPa) 

85.0 

Damage  constants 

D i 

0.04 

d2 

1.0 

Efmmh 

0.01 

Pressure  constants 

^crush  (GPa) 

0.016 

t1  crush 

0.001 

Aock  (GPa) 

0.80 

A^lock 

0.10 

K\  (GPa) 

85.0 

K2  (GPa) 

-171.0 

K3  (GPa) 

208.0 

T  (GPa) 

0.004 

Strength  constants 

A 

0.79 

B 

1.60 

C 

0.007 

N 

0.61 

/'  (GPa) 

0.048 

5max 

7.0 

k  (s-1) 

10“6 

Unconfined  compressive  strength  of  concrete  =  48  MPa. 
Ermin  is  the  amount  of  plastic  strain  before  fracture. 


Table  1.  For  comparisons,  the  computations  are  also  car¬ 
ried  out  for  the  L-E  model.  The  variations  of  peak  axial 
stress  oy  with  relative  distance  y  (refer  to  the  definition 
shown  in  Fig.  1)  are  shown  in  Fig.  4(a)  for  the  J-H-C 
model  and  Fig.  4(b)  for  the  L-E  model.  In  both  cases, 
the  peak  oy  decreases  with  increase  in  radius  R  for  a  spec¬ 
ified  position  y.  On  the  other  hand,  the  peak  oy  increases 
with  relative  distance  y  for  a  nonzero  R.  Comparing  the 
magnitudes  of  the  peak  axial  stress  oy  in  Fig.  4(a)  and 
(b),  it  can  be  seen  that  the  elastoplastic  properties  of  con¬ 
crete  plays  a  crucial  role  on  the  magnitude  of  peak  axial 
stress.  For  instance,  when  R  =  0.30  m  and  y  =  0.60  m, 
the  peak  oy  is  51.8  MPa  for  the  L-E  model  and  only 
25.2  MPa  for  the  J-H-C  model.  The  change  of  decay  factor 
DF  with  R  is  shown  in  Fig.  5  for  the  J-H-C  model.  This 
figure  shows  that  the  DF  increases  with  R  but  reduces  with 
increasing  y. 

The  hydrostatic  pressure  contours  are  shown  in 
Fig.  6(a)  and  the  damage  contours  are  shown  in 
Fig.  6(b).  These  figures  show  that  the  stress-wave  is  dis¬ 
rupted  due  to  the  presence  of  cavity.  A  stress  concentra¬ 
tion  is  observed  in  the  left  and  right  sides  of  the  cavity, 
with  a  maximum  value  of  almost  12  MPa.  A  zero-stress 
region  is  noted  in  the  adjacent  domain  beneath  the  cavity. 
The  peak  oy  beneath  the  cavity  is  evidently  lower.  Refer¬ 
ring  to  Fig.  6  (b),  the  damage  zone  is  mainly  located  at 
the  upper  surface  of  the  problem  domain.  Therefore,  the 
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Fig.  4.  Effect  of  different  radii  R  of  circular  cavity  on  peak  axial  stress  oy. 
(a)  J-H-C  model,  (b)  L-E  model. 


Fig.  5.  Variation  of  decay  factors  for  the  peak  oy  with  different  radii  R  of 
cavity  for  JHC  model. 


air-blast  loading  may  cause  spalling  at  the  roof  of  the  cav¬ 
ity.  Special  protection  is  necessary  when  the  peak  value  of 
blast  loading  is  much  larger. 


4.1.2.  Effect  of  cavity  position  e 

The  effect  of  the  cavity  position  e  on  the  attenuation  of 
stress  waves  is  shown  in  Fig.  7.  When  R  is  fixed,  the  decay 
factors  DF  are  practically  identical  for  all  e  values.  This 
phenomenon  can  be  also  observed  from  the  time-hydro- 
static  pressure  curves  at  y  =  0.5  m  shown  in  Fig.  8.  The 
peak  value  of  hydrostatic  pressure  does  not  vary  much 
for  all  e  values.  It  should  be  noted  that  the  above  results 
are  obtained  when  the  cavity  has  not  collapsed. 

4.1.3.  Empirical  formula  for  decay  factor 

Based  on  the  above  numerical  results,  it  is  established 
that  DF  is  essentially  a  function  of  R  and  y: 

DF  =f(R,  y)  (9) 

where  the  relative  distance  y  =  Y  —  2R  —  e. 

In  order  to  determine  the  particular  form  of  Eq.  (9), 
one  value  of  e  (0.30  m),  four  values  of  R  (0.15,  0.30, 
0.45  and  0.60  m)  and  five  values  of  y  (0.2,  0.4,  0.6,  0.8 
and  1.0  m)  are  chosen  from  Fig.  5.  Through  a  curve  fit¬ 
ting  regression  analysis,  the  DF  for  the  peak  oy  at  a  spec¬ 
ified  point  on  the  F-axis  is  approximately  expressed  as 
follows: 

DF  =  1  -  exp  (-  (10) 

This  formula  shows  that  DF  increases  with  R  and 
decreases  with  increasing  y.  When  |  approaches  to  zero, 
the  DF  is  almost  equal  to  unity.  Fig.  9  shows  the  compar¬ 
ison  of  numerical  results  with  the  predictions  from  the 
above  empirical  formula.  It  is  clear  that  this  empirical  for¬ 
mula  gives  a  reasonable  good  fitting  to  the  numerical 
results.  Therefore,  Eq.  (10)  can  be  used  to  estimate  the 
radius  R  of  a  circular  cavity  if  the  expected  DF  (at  a  rela¬ 
tive  distance  y)  and  position  e  are  known.  On  the  other 
hand,  when  DF  and  R  are  given,  the  relative  distance  y 
can  be  also  determined  from  Eq.  (10). 

4.2.  Rectangular  cavity 

Rectangular  cavities  are  also  often  built  in  the  defense 
layer.  The  effect  of  a  rectangular  cavity  on  the  attenuation 
of  stress-waves  will  be  numerically  investigated  in  this  sec¬ 
tion.  The  embedded  cavity  has  the  dimensions  of  2/  in 
length  and  h  in  width.  Other  dimensions  are  indicated 
in  Fig.  10.  The  dimension  Y  has  the  following 
relationship: 

Y  =y  +  e  +  h  (11) 

In  the  computations,  the  constitutive  model  of  concrete 
and  boundary  conditions  are  the  same  as  those  for  the  cir¬ 
cular  cavity  case.  A  half  domain  is  also  considered  here 
owing  to  symmetry. 

4.2.1.  Effect  of  cavity  length  l 

Five  lengths  are  taken  in  the  computations:  0.0,  0.3,  0.6, 
0.9  and  1.2  m.  Other  geometrical  parameters  are  fixed  as 


Z.-l.  Wang  et  al.  /  Computers  and  Geotechnics  33  (2006)  132-141 


137 


Time  =  439.97 


Fringe  Levels 
1.921e-04 
1.689e-04  I 
1.457e-04  _l 
1.224e-04  _ 
9.924e  05  _ 
7.603e-05  _ 
5.283e-05  . 
2.962e-05  _ 
6.41 4e-06 
-1.679e-05  _ 
-4.000e-05  I 


a 


Time  =  439.97 


b 


Fringe  Levels 
4.093e-02 
3.683e-02 
3.274e-02 
2.865e*02 
2.456e-02 
2.046e-02 
1.637e-02 
1.228e-02 
8.185e-03 
4.093e  03 
0.000e*00 


Fig.  6.  Contours  of  pressure  and  damage  for  R  =  0.30  m  cavity  at  t  =  440  ps  when  subject  to  air-blast  wave,  (a)  Pressure  (unit:  Mbar.  1  Mbar  =  105  MPa), 
(b)  Damage. 


Fig.  7.  Variation  of  decay  factors  for  peak  ov  with  different  positions  e  of 
cavity. 


e  =  0.6m  and  h  =  0.3  m.  Fig.  1 1  shows  the  effect  of  differ¬ 
ent  lengths  /  on  the  attenuation  of  the  peak  oy.  The  results 
show  that  the  peak  oy  increases  with  relative  distance  y  and 
decreases  with  increasing  length  /.  The  corresponding  rela¬ 
tionship  between  DF  for  the  peak  oy  and  relative  distance  y 
is  shown  in  Fig.  12.  It  can  be  established  that  DF  decreases 
with  increasing  y  and  increases  with  cavity  length  for  the 
same  y. 


Fig.  8.  Effect  of  cavity  position  e  on  hydrostatic  pressure  history  at 
y  =  0.5  m  point  on  7-axis. 


4.2.2.  Effect  of  cavity  width  h 

Four  values  of  h  =  0.3,  0.6,  0.9,  1.2  m  when 

l=e  =  0.3  m  are  considered  to  study  the  effect  of  cavity 
width  on  the  attenuation  of  stress- waves.  The  numerical 
results  are  shown  in  Fig.  13.  It  is  noted  that  the  cavity 
width  h  has  little  influence  on  the  decay  factor  for  the 
peak  Gy.  Fig.  14  shows  the  damage  contours  for 
h  =  0.60  m  at  a  time  of  524  ps,  which  indicates  that 
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Fig.  9.  Comparison  of  empirical  formulae  and  numerical  results  of  decay 
factors  for  the  peak  oy. 


Fig.  12.  Variation  of  decay  factors  for  the  peak  oy  with  different  lengths  / 
of  cavity. 


Air-blast  pressure  Q 


Fig.  10.  Computational  domain  and  dimensions  for  rectangular  cavity 
scheme. 


Fig.  13.  Variation  of  decay  factors  for  peak  oy  with  different  widths  h  of 
cavity. 
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Fig.  11.  Variation  of  the  peak  oy  with  different  lengths  /  of  cavity. 

the  high  damage  zone  is  located  only  near  the  upper 
surface  and  approaches  zero  at  the  zone  beneath  the 
cavity. 


4.2.3.  Effect  of  cavity  position  e 

The  effect  of  cavity  position  e  is  examined  here.  The  e 
values  are  taken  as  0.3,  0.6,  0.9  and  1.2  m  and 
1  =  h  =  0.3  m.  Fig.  15  shows  the  decay  factor  of  the  peak 
oy  versus  relative  distance  y.  It  is  evident  that  all  the  curves 
are  practically  identical.  Therefore,  the  position  e  of  the 
cavity  has  little  influence  on  the  decay  factor  for  the  peak 
Oy,  similar  to  the  case  of  circular  cavity  presented  earlier. 
Fig.  16  compares  the  curves  of  oy  versus  t.  For  all  e  values, 
the  peak  values  do  not  vary  much,  but  the  arrival  time  of 
the  peak  values  gradually  lags  with  an  increase  of  the  prop¬ 
agation  distance  of  stress- waves. 

4.2.4.  Empirical  formula  for  decay  factor 

Based  on  the  above  numerical  results,  the  effects  of  h  and 
e  on  the  peak  oy  are  established  to  be  insignificant.  Thus, 
DF  is  only  a  function  of  length  /  and  relative  distance  y 

DF  =f(l,y)  (12) 
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Fig.  14.  Damage  distribution  for  the  h  =  0.60  m  (e  =  /  =  0.30  m)  case  at  t  =  524  jis. 


Fig.  15.  Variation  of  the  decay  factors  for  peak  oy  with  different  positions 
e  of  cavity. 


Time  (pis) 


Fig.  16.  Effect  of  cavity  position  e  on  oy  history  at  y  =  0.8  m  point  on  Y- 
axis. 


Based  on  the  results  shown  in  Fig.  12,  one  value  of  h 
(0.3  m),  four  values  of  /  (0.3,  0.6,  0.9  and  1.2  m)  and  five 
values  of  y  (0.2,  0.4,  0.6,  0.8  and  1.0  m)  are  chosen  in  the 
computations  to  obtain  the  following  regression  curve  for 
the  peak  oy  at  the  points  just  beneath  the  cavity 


DF  =  exp 


(13) 


This  empirical  formula  shows  that  DF  increases  with  / 
but  decreases  with  increasing  y.  DF  is  equal  to  unity  when 
y  approaches  0.  Fig.  17  shows  a  comparison  between  the 
numerical  results  and  the  predictions  from  the  empirical 
formula.  This  empirical  formula  is  in  good  agreement  with 
the  numerical  results.  Therefore,  Eq.  (13)  can  be  used  to 
estimate  the  length  /  (or  the  whole  length  21)  of  a  rectangu¬ 
lar  cavity  if  the  expected  DF  at  a  relative  distance  y  is 
given.  When  DF  and  /  are  given,  Eq.  (13)  can  be  also  used 
to  estimate  the  nearest  safe  distance  (i.e.,  ymin). 

In  defense  engineering,  the  DF  value  at  a  specified  posi¬ 
tion  beneath  a  cavity  is  required  to  ensure  the  safety  of  the 
underground  structures.  For  this  purpose,  Eqs.  (10)  and 
(13)  can  be  used  to  estimate  the  minimum  radius  (Rm in) 
or  length  (/min)  of  a  cavity.  For  example,  when  y  =  1.0  m 
and  DF  ^  0.5,  /min  should  be  at  least  0.40  m  for  a  rectan¬ 
gular  cavity  while  the  radius  7?min  should  be  at  least  0.69  m 
for  a  circular  cavity.  The  attenuation  effect  of  a  circular 
cavity  is  much  lower  than  that  of  a  rectangular  cavity  if 
R  =  l. 


Fig.  17.  Comparison  of  empirical  formulae  and  numerical  results  of  the 
decay  factors  for  the  peak  oy. 
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Fig.  18.  Computational  domain  for  multiple  rectangular  cavities  scheme  (unit:  cm). 


4.3.  Effects  of  multiple  cavities 

The  above  results  indicate  that  a  single  cavity  exerts  a 
great  influence  on  the  attenuation  of  stress- waves.  In  prac¬ 
tice,  a  screen  of  multiple  artificial  cavities  is  usually  formed 
to  insulate  and  diffract  stress-waves.  Fig.  18  depicts  a  typ- 


Fig.  19.  Comparisons  of  screen  effect  between  single  and  multiple 
rectangular  cavities,  (a)  Curves  of  axial  stress  versus  relative  distance, 
(b)  Curves  of  axial  stress  versus  lateral  distance. 


ical  defense  layer  embedded  with  five  rectangular  cavities 
(0.9  x  0.3  m).  The  spacing  between  two  neighborhood  cav¬ 
ities  for  this  case  study  is  0.25  m.  The  same  triangular 
impulse  loading  shown  in  Fig.  2  is  applied  to  the  upper  sur¬ 
face.  The  insulating  effect  between  multiple  cavities  and  a 
single  cavity  on  the  reduction  of  the  peak  oy  when 
e  =  0.9  m  and  h  —  0.3  m  is  compared.  Fig.  19(a)  shows  that 
the  decay  factor  for  the  peak  oy  of  multiple  cavities  is  much 
higher  than  that  of  single  cavity.  However,  when  the  rela¬ 
tive  distance  is  less  than  1.0  m,  the  effect  of  multiple  cavities 
(/  =  0.45  m)  is  rather  close  to  that  of  single  cavity  with 
/=  0.6  m.  Fig.  19(b)  compares  the  decay  factor  for  the 
peak  oy  in  the  lateral  direction  at  the  points  for 
y  =  1.0  m.  With  the  increase  of  lateral  distance  x,  the  peak 
oy  of  multiple  cavities  fluctuates  around  a  stable  value  (i.e., 
12.0  MPa).  The  peak  value  is  also  much  lower  than  that  of 
single  cavity.  Besides,  one  can  see  that  the  peak  oy  of  single 
cavity  is  almost  the  same  to  the  one  without  any  cavity  (i.e., 
38.2  MPa)  when  x  >  1.2  m.  Apparently,  the  shielding  effect 
is  much  better  for  multiple  cavities  than  for  a  single  cavity. 
A  screen  against  stress-waves  is  formed  through  multiple 
cavities.  This  screen  can  effectively  reduce  the  magnitude 
of  stress- waves  beneath  the  cavity. 


5.  Conclusions 

The  effects  of  circular  and  rectangular  cavities  in  a  con¬ 
crete  defense  layer  on  the  diffraction  and  attenuation  of 
stress-waves  are  investigated  using  the  J-H-C  model. 
Two  empirical  formulae  are  proposed  to  describe  the  decay 
factor  of  the  peak  axial  stress.  In  addition,  the  screen  effect 
due  to  multiple  rectangular  cavities  is  analyzed.  From  these 
studies,  the  following  conclusions  and  understandings  can 
be  drawn. 

(1)  The  elastoplastic  and  damage  properties  of  concrete 
have  great  influence  on  the  magnitude  of  axial  stress 
oy  in  the  concrete  defense  layer.  The  J-H-C  model  is 
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capable  of  evaluating  the  elastoplastic  and  damage 
effects  of  concrete  on  the  reduction  of  the  axial  stress 
Gy.  Thus,  the  model  can  be  employed  for  the  design  of 
concrete  defense  layer. 

(2)  Cavities  embedded  in  the  concrete  defense  layer  have 
a  vital  effect  on  the  attenuation  of  stress- waves.  For  a 
circular  cavity,  the  radius  R  has  a  significant  effect  on 
the  reduction  of  stress- waves.  For  a  rectangular  cav¬ 
ity,  the  decay  factor  largely  depends  on  its  length, 
while  its  width  has  little  effect  if  the  width  is  larger 
than  0.3  m  in  our  computations. 

(3)  Empirical  formulae  are  obtained  to  correlate  the 
decay  factor,  dimensions  of  a  cavity  and  relative  dis¬ 
tance.  These  formulae  can  be  employed  to  determine 
the  minimum  length  or  radius  of  a  cavity,  or  mini¬ 
mum  safe  distance  ymin  when  the  decay  factor  is 
given.  Therefore,  these  formulae  can  serve  as  useful 
tools  in  the  design  of  intelligent  defense  layer  in 
defense  engineering. 

(4)  Multiple  cavities  can  form  a  screen  to  reduce  the 
magnitude  of  stress- waves.  The  decay  factor  at  any 
point  beneath  this  screen  is  much  higher  than  that 
of  a  single  cavity.  The  shielding  effect  is  also  pro¬ 
nounced  especially  in  the  transverse  direction. 

It  should  be  noted  that  the  above  conclusions  are  appli¬ 
cable  only  when  the  damage  of  concrete  does  not  initiate 
the  collapse  of  the  cavity.  Further  work  is  necessary  when 
the  damage  of  concrete  defense  layer  is  severe  enough  to 
cause  a  complete  collapse  of  the  cavity. 
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